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Abstract 



In this article, we consider abstract linear conservative systems and 
their time-discrete counterparts. Our main result is a representation for- 
mula expressing solutions of the continuous model through the solution 

?«>^ I of the corresponding time-discrete one. As an application, we show how 

observability properties for the time continuous model yield uniform (with 
respect to the time-step) observability results for its time-discrete approx- 
imation counterparts, provided the initial data are suitably filtered. The 
main output of this approach is the estimate on the time under which we 
can guarantee uniform observability for the time-discrete models. Besides, 
using a reverse representation formula, we also prove that this estimate 

k>( ' on the time of uniform observability for the time-discrete models is sharp. 

'%_] ' We then conclude with some general comments and open problems. 
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1 Introduction 

1.1 Setting 



Assume that A is a skew-adjoint unbounded operator defined in a Hilbert space 
X with dense domain T>{A) C X and compact resolvent. We consider the 
equation 

y'^Ay, teR, y{Q) = y°^X. (1.1) 

This equation is well-posed for t G R and, A being skew-symmetric, solutions 
y{t) of (|l.ip have constant norms ||j/(i)||_:t.. Often ||y(i)||_;^, /2 is referred to as 
the energy of the system. 

Several classical systems fit this abstract setting, as for instance the wave 
equation, Schrodinger's and Maxwell's equations, among many others. 

Our primary goal is to establish a link between the solutions y of (jl.l[) and 
the solutions of time-discrete versions of (|l.ip . Since the solutions y of (II. ip 
have constant energy, we will reduce our analysis to time-discretization schemes 
which preserve the A'-norm \\-\\x of the solutions. 

To state our results precisely, we need to introduce filtered spaces of solu- 
tions. For, we write down the spectral decomposition of the operator A, which 
is given by a sequence of purely imaginary eigenvalues («/Xj), corresponding to 
an orthonormal basis of X constituted by the eigenvectors $j. For (5 > 0, we 
define the filtered space 

<Z{5) = span{$j corresponding to eigenvalues |^j| < 5}. (1.2) 

One easily checks that this space is left invariant by the equation (|l.ip , and thus 
we will identify the space of trajectories y solutions of (jl.ip lying in £((5) with 
the space of initial data y° lying in €.{5). 

Let us now describe the time-discretization schemes under consideration. 
We assume that the time-discretization of system (|l.ip discretized with a 
time-step r > takes the form 

y^+^^Try";, fceZ, yl = y'. (1.3) 

Here, y^ denotes an approximation of the solution y of (jl.ip at time fcr, and T,- 
is assumed to be an approximation of cxp(T^). 

To be more precise, we assume that T,- is a linear operator which has the 
same eigenvectors as the operator A, and such that, for some function /, 

T,-$j = cxpiifij^rT)^]., where jjlj^t ~ —/{^jt). (1-4) 

r 

Actually, this operator may also be written as 

T, = exp (^/(-^^T)) . 

Let us now make precise the assumptions we impose on /. 
First, we assume that / is C°° and satisfies 

/(0)-0, /'(0) = 1. (1.5) 



This assumption on / is satisfied for most of the time-discretization approxi- 
mation schemes: Roughly speaking, this is equivalent to the consistency of the 
time-discretization. 

To simplify notations and avoid technical developments, we furthermore as- 
sume that / is odd. This is not necessary in our arguments, but in practice, 
most of the time-discretization schemes modeled by (I1.3P fit in this class. 

We also assume that / is real valued and 

f:{~R,R)^{-7T,7r), (1.6) 

where R E M!i_U{oo} may be infinite. The fact that / is real- valued is equivalent 
to the fact that the norms of solutions of (jl.3p are constant in the discrete time 
fc G Z and ensures the stability of the numerical scheme. The range of / is 
limited to (— 7r,7r) to avoid aliasing since one cannot measure oscillations at 
frequencies higher than vr/r on a discrete mesh of mesh-size r. 

In particular, for (11.31) to be well-defined, we always consider solutions of 
(II. 3p with initial data lying in some class of filtered data at the scale r, namely 
(L{S/t) for some 5 g (0,i?). 

We also assume that for all S < R, 

inf{/'(a); \a\ < S} > (1.7) 

to ensure the invertibility of /. 

We finally define g : {—f{R),f{R)) — > M as the inverse function of / on 
{—R,R). This can be done due to assumption (|1.7p . 

Before going further, let us point out that several classical time-discretization 
schemes fit the abstract setting provided by Assumptions (|1.5p - (|1.6p - (|1.7p . see 
Section [5] for several examples. 



1.2 A representation formula 

We are now in position to state the following result, whose proof can be found 
in Section [21 

Theorem 1.1. Let f be a smooth function describing the time discrete operator 
Tr as in pTil) . assume ([T3)) ^ PT^ - PTT)) . and fix S e {0,R). 

Let X be a C°° function compactly supported in (— /(i?), /(i?)) and equal to 
1 in (-/(<5), f{5)). For r > 0, define then pr{t, s)fort£R, s e K as 

Pr{t, s)^:^ r'^ exp f-M^^) ^ (^^r) e^^^^ d^.- (1-8) 



2^ J-./r 

Then, if ip G £((5/r) and y^ is the corresponding solution of (jl.3p . the function 
y{t) defined by 

y{t)^TY,Pr{t,kT)y^r (1-9) 

fcez 
is the solution of (|l.ip with initial data y^ . 

Remark 1.2. Note that when f is bijective from {—R,R) to (— 7r,7r), the above 
representation formula does not require the use of the cut-off function x and 
allows describing all solutions y of (jl.ip for all initial data in y'^ e (^{R/t), see 
Remark \3.S\ for further details. 



We now give some informal arguments motivating the representation formula 
in Theorem II .11 

First, we remark that under the assumptions on the operator A, solutions y 
of (jl.lD admit the following spectral decomposition: 

yit) = J2a,e^^^%, (1.10) 

j 

where (oj ) are the coefRcients of the initial datum on the basis $j : 

j 

Similarly, due to (11.41) , the solutions yr of the time-discrete system (|1.3p with 
initial datum y'^ as in (|l.lip can be written as 

2/^- = Ea,e^^-'=-a>„ (1.12) 

j 

where the iij^r are defined by p.4p . Obviously, we can then extrapolate a time- 
continuous version Zr of yr : 

z.(s) = Ea,e^''-^<I>,. (1.13) 

j 

Here, we use the notation s for the continuous time-variable corresponding to 
the a priori time-discrete dynamics, and z^ for the new state variable to avoid 
confusion with time t and state y of the continuous dynamics. 

Transmutation or subordination refers to the possibility of expressing one 
semi-group as a function of another one. In our context, this simply consists 
in writing the solutions y of (jl.ip as a function of Zt(s) by means of a kernel 
p — p{t, s) under the form 

y{t)= Pr{t,s)Zris)ds, t £ R. (1.14) 

For this to be done, using the explicit expressions (ll.lOp and (|1.13p . the kernel 
p{t, s) has to be built so that 



gip-jt _ 



I Pr{t,s)e'''^--'ds, teR,peR. (1.15) 

Jr 



Interpreting the right-hand side of p.lSp as a Fourier transform in s and taking 
into account that iij^r — fifJ-j''')/''' with / invertible, naturally leads to 

^^^^''^ ^ ^I'^'P (^„!5(^^ e^..rf^. (1.16) 

This formula is a simplified version of the one of p-r in (jl.Sp . in which the above 
time-continuous function z^ in (|1.13p has to be replaced by the time-discrete 
function y^ in (|1.12p and the filtering operator x has been introduced. 



1.3 Application to observability 



As an application of Theorem 11.11 we consider an observation problem for the 
equation (ll.ip and its time-discrete counterparts for (ll.3p . This problem, in- 
spired in control theoretical issues, concerns the possibility of recovering the 
full energy of solutions out of partial measurements. The question is relevant 
both in the continuous and in the time-discrete frame, and in the later, a natu- 
ral question arising in numerical analysis is to know whether the observability 
property is uniform with respect to the time-step. Indeed, this problem is the 
dual version of the classical controllability one and the uniformity (with respect 
to time-step) of the observability inequality is equivalent to the convergence of 
ntimerical controls towards the continuous ones as the time-step tends to zero 
(see [501 [IS]). 

We thus consider an observation operator B taking value in some Hilbert 
space 14 and assumed to belong to £,{'D{Ap),U) for some p £ N. To be more 
precise, we assume that there exists Cp > such that 

\\By\\u < Cl {\\A^y\\l + ||y||^) , Vy e ViA^. (1.17) 

We furthermore assume that equation (II. ip is observable through B at time 
To > 0, meaning that there exists a constant Co such that, for all y" e 25(^1^), 
the solution y{t) of (II. ip with initial data y° satisfies 

\\y'f^<C,r\\By{t)\\l,dt. (1.18) 

Estimate (|1.18p is the so-called observability estimate for (II. ip . As it has been 
established in the classical works [3 [IS] j this property is essentially equivalent 
to the controllability of the adjoint system: 

y', = Avc + B*v, ie(0,r), yM = vl (1.19) 

where t; is a control function in L^{0,T;U). Actually the control for (J1.19p 
can be built by minimizing a suitable quadratic functional over the class of 
solutions (jl.ip . We do not give details on these links in this article and we refer 
the interested reader to (23] [35] ■ Also note that the observability property is 
also closely linked to some inverse problems, see e.g. the work W. for precise 
statements in a setting similar to ours and the references therein. 

We are thus interested in deriving discrete versions of (|1.18p for solutions y^. 
of (|1.3I) . Namely, we are asking if, given 6 G (0, i?), there exist a time T > and 
constants C > and To > independent of r > such that for all t e (0,to), 
solutions j/t- of (|1.3p with initial data y° S £((5/r) satisfy 

\\yll<(^^ E WrWl- (1-20) 

fcTe(o,T) 

This is a discrete observability estimate for the time-discrete equations (|1.3p 
which is uniform with respect to the time-discretization parameter, provided the 
initial state is suitably filtered. As it is well-known, see e.g. the survey articles 
[501 [T5] , this is needed to derive algorithms to compute convergent sequences of 
discrete controls approximating the control of the continuous dynamics. 

Our second main result, whose proof is given in Section [3J is the following 



Theorem 1.3. Assume that B G £(I'(.4p),W) for some p E N andU an Hilhert 
space, and that B satisfies (|1.17p with constant Cp. 

Assume that equation (jl.ip is observable through B at time Tq with constant 
Co, i.e. for all y^ G ^[A^), the solution y[t) of (|l.ip with initial data y" 
satisfies (jl.181) . 

Let f be a smooth function describing the time discrete operator Tr as in 
([Li)) , assume p3]) - pT^ -PT7 )) . and fix 5 G (0,i?). Then, for all 

there exist positive constants C and tq > depending on f, T, 5, p, Cp, Cq, Tq such 
that for any r G (0,tq), solutions yr of (jl.3p lying in £((5/t) satisfy (|1.18|) . 



Theorem 11.31 was already proved in [13] under an additional admissibility 
property for the operator B (see [5S] for a definition, or Section 14.21 in our 
context), but with an estimate on the observability time which is worse than the 
one we prove here in (jl.2ip . Indeed, the technique in [T3], based on a resolvent 
characterization of observability due to [H I33| , does not yield explicit estimates 
on the time of observability. A variant of this strategy was developed in [TTI for 
time-discrete approximations of Schrodinger equations in a geometric setting 
in which the corresponding wave equation is observable, where it was proved 
that the time-discrete approximations of Schrodinger equations are uniformly 
observable in any time T > 0. In that case indeed, the resolvent estimates 
behave much better at high-frequencies, see e.g. 49], and the low frequency 
components of the solutions can be handled following the arguments of |18) . 
These resolvent estimates were also used to derive observability estimates for 
space semi-discrete conservative systems in filtered classes of initial data, see 
pm [T^ ;36: , but still with non-explicit estimates on the time of observability. 

The new ingredient introduced in this article that allows us to improve the 
results in |13j is the representation formula (jl.9p and careful estimates on the 
kernel pr in (jl.Sp . This will be done using classical techniques of harmonic 
analysis, and in particular the oscillatory phase lemma. In particular. Proposi- 
tion [2il] shows that, for all e > 0, pr is polynomially small with respect to r > 
to any arbitrary order in the set 

|(t,s) G (0,T) xR, s.t. t + e < s inf {/'(a)} or s sup {f'{a)} < t - e]. 

I \a\<S+e \a\<S+e -' 

To give some insights on this result, let us again consider the informal arguments 
given above and remark that Zt{s) in (11.131) formally solves 

where 9J is the operator defined on the Fourier basis by 

V^iGM, 5J (e^''") = i^^e'^^ (1.22) 

Thus, the kernel pr in p.l6p formally satisfies the transport-like equation 

dtPr - dlpr = 0, (t, s) G K X M, (1.23) 



with an initial data pr{t = 0, s) = 5q{s) where 5q{s) is the Dirac function, since, 
according to (|1.16p . pr{t = 0, •) is the Fourier transform of the function taking 
value 1 identically. 

Note that the solution of (9t — 9s)p* = with initial data p*(t = 0) = So{s) is 
simply the Dirac delta transported at velocity one: /5*(i, s) — St{s). Of course, 
for the kernel under consideration, 9J is not a classical differential operator, 
but it is nevertheless very close to the classical derivation operator dg at low 
frequency since g'{0) = 1. We can thus interpret (|1.23p as a transport equation 
with some added dispersion effects at frequencies of the order of l/r. 

Let us also point out that the representation formula in (jl.9p given through 
the kernel p^ in (jl.Sp is not far from a Fourier integral operator - see e.g. 
P5] - with the phase ip{t,s,pr) — fJ-rS — g{fiT)t, and similar techniques can 
be employed. In particular the above localization result can also be seen as a 
counterpart of the fact that the kernel of a Fourier integral operator with phase 
ip is regularizing outside of the set in which dfj_^p — 0. 

Next, in Theorem 14. 11 we discuss a reverse representation formula giving so- 
lutions of the time-discrete equation (|1.3p in terms of solutions of the continuous 
equation (|l.ip . This allows us to prove admissibility results for solutions of (II. 3p 
uniformly with respect to the time-discretization parameter r, see Theorem 14. 5 [ 
and the optimality of the time estimates in (|1.2ip , see Section 14.31 

Let us also point out that, similarly as in [13], our approach can also be 
applied in the context of fully-discrete schemes: Indeed, to derive observability 
estimates for fully discrete approximations of (II. ip that are uniform in both time 
and space discretization parameters, our approach shows that it is sufficient 
to prove observability estimates for time-continuous and space semi-discrete 
approximation schemes that are uniform in the space discretization parameter, 
see Section 15.11 for precise statements. Also note that our results can be used 
to recover discrete Ingham inequalities in a slightly different setting as the one 
in [JU], see Section [^21 We also explain how our strategy applies in the case of 
weak observability estimates in Section [ 



1.4 Related results 

Our approach is inspired by several previous works which establish represen- 
tation formula for solutions of one equation through the solution of another 
equation. As we have already said, this technique is called subordination (in 
particular in the context of functional analysis, see e.g. [13] and references 
therein) or transmutation. For instance, the work [21 provides estimates on 
the heat kernel thanks to an analysis of the corresponding wave equation and 
the so-called Kannai transform expressing solutions of the heat equation in 
terms of solutions of the wave equation, and the work [22^ proposes a study of 
singular problems based on a representation formula adding one dimension to 
the problem. 

In the context of control theory in which we focus here, the so-called Fourier 
Bros lagnoniltzer transform, linking solutions of the wave equation to a suitable 
elliptic operator, provides an efficient tool to prove quantified unique continua- 
tion properties, see e.g. [1511321132] ■ 

Similarly, using a suitable transformation linking the wave and heat equa- 
tions. Miller in [331 ES] derived estimates on the cost of controllability of the 



heat equations in small time, later improved in |47| . Similar estimates, related 
to the characterization of the reachable set for heat equations, can be found in 
[H] , The main difference between the transforms in [Ml [3S] and [Tl] is that, 
whereas the articles [3H ES] express solutions of the heat equations in terms of 
solutions of the wave equation, [M] is based on a transform expressing solutions 
of the wave equation in terms of the solutions of the heat equation. Note that 
the robustness of the transformation of [T3] is illustrated by the fact that the 
weak observation properties derived in |42) for the waves in general geometric 
setting (not necessarily satisfying the geometric control conditions) can be used 
to recover control properties for the heat operator. Let us also point out that 
these transmutation techniques can also be used to derive numerical schemes to 
compute approximate controls for the heat equations |37] . 

We also emphasize that Theorem ll.3l states observability results for the time- 
discrete schemes (jl.41) observed through an observation operator B provided the 
continuous model (jl.ll) is observable through B. Hence the first step is to check 
the observability property for the continuous model, and these observability 
properties have to be checked in each situation. When the continuous model 
(jl.ip stands for a wave equation and 6 is a distributed observation operator or 
a boundary observation, the necessary and sufficient condition for observabil- 
ity is the so-called geometric control condition, see [21 13] • For what concerns 
plate or Schrodinger's equations, the geometry still plays a role, but due to the 
infinite speed of propagation the situation is more intricate: for these models, 
observability holds when the geometric control condition is satisfied (see e.g. 
[22] where this result is derived using transmutation techniques) , but other less 
restrictive geometric settings still enjoy observability properties, see e.g. [1]. 

1.5 Outline 

The article is organized as follows. In Section |21 we give several instances of 
classical time-discretization schemes that fit the abstract setting (jl.4p with a 
function / satisfying Assumptions (|1.5p - (|1.6p - (|1.7p . Section [31 is devoted to the 
proofs of Theorem 11.11 and Theorem II. 31 In Section |3] we present a reverse rep- 
resentation formula and discuss its application to uniform admissibility results 
and use it to prove the optimality of the time estimate (ll.2ip . In Section jS] we 
present some further comments. In Section [6l we end up discussing some open 
problems. 

2 Some admissible time-discretization schemes 

In this section, we provide several classical time-discretization schemes that fit 
the setting (fO|) - (fri|) and satisfy assumptions (|T3)) - prB)) - pr7|) . 

2.1 The midpoint scheme 

Perhaps the simplest time-discretization of (|l.ip which preserves the energy is 
the midpoint scheme: for t > 0, the time-discrete equation is given by: 

y^^^^jl^^/y^+_^\ fceZ, yO-y". (2.1) 



Thus, if ^j is an eigenvector of A with eigenvalue i/ij, the solution yr of 
(|2.ip with initial data y^ — $j is given by 

where fj,j,r is defined by exp{i^j^rT) = r-^ — j^, 

1 — lfljT/2 

yielding 

fj.j_TT — 2 arctan(/ijr/2). 

Hence the midpoint scheme (j2.ip fits the assumptions of Theorem 11.31 by 
setting /(a) ~ 2arctan(a/2) and thus R = 00. 

2.2 The fourth order Gauss Method 

Let us present the so-called fourth order Gauss method for discretizing (jl.lL 
which enters the frame of Runge-Kutta methods, see for instance [T7] . 
It reads as follows: 



^ 1 -^^ 



Kf =^ U^+T^ayKM , i = l,2 

iA'+l - I/'-' -I- -I'f,-'' -t- k'^'I / 1 1 v. 

Vt -Vt + -^\'^1+>^2), (ry.\-\ 4 4^-6 

. Vr^y ^^ given, V 4 + 6 4 



(2.2) 



An easy computation shows that, for any A; G Z, 

\{4 + 4)^[ld-^ + ^-^j Ay";, 
which allows to rewrite 

The spectral decomposition of the semi-discrete scheme (j2.2p can easily be 
performed. If y^ = '^j, we obtain 

y^ = exp I i — fcr ) y°, where /(a) = 2 arctan I 

\ T / \ 12 — a-^ 

Hence the discretization (|2.2p fits the assumption of Theorem 11.31 by setting 

/ : (-2%/3, 2%/3) ^ M; /(a) = 2arctan ( — ^ , 

\ 12 — a^ 

Note that, here, the function / is limited to the range where R = 2\/3- 



2.3 The Newmark method for second order in time equa- 
tions 

The Newmark method is designed for second order in time equations such as, in 
particular, the wave equation. Namely, let ^o be a self-adjoint operator defined 
on an Hilbert space Xq with dense domain ^{Ao) and with compact resolvent, 
and consider the following equation: 

ip"+Ao^ = 0, teR, {ip,ip'm^i^°,ip^)eViAl^^)xXo. (2.3) 

The Newmark method yields, for j3 e [0, 1/4], the following time-discrete scheme: 

+ ^0 (/3^^+i + (1 - 2(3)^'; + 13^';-') = 0, 



^k+1 



^^-1-2^^ 



V°r+Vl ^l~V>r 



(/,/)eA'2. 



2 ' r 
System (|2.4p is conservative and preserves the discrete energy: 



(2.4) 



^fc+l/2 ^ 



^1/2 f ^'r + Vr^' 



LP'; 



fe+i 



Vr 



Xo 



+ (4/3-1)- 



Xq 



.1/2 /'(^^+^-(^^ 



System fO)) fits the abstract setting pTTI) by setting X = 'D{Aq' ) x Xq, 



Xo 
l/2x 



(2.5) 



i1/2n 



,1/2 ,1/2 



where V^Aj ) is endowed with the scalar product {Aq -^Aq ■)xa £^nd 

A- 



/ 

-Ao 



But another way to write (j2.3p under the form (jl.ip is to set A" = A^q , 



a/2 



yi = </? + z^q' (p, y2 = if - iAq If 



1/2, 



y^ = 



iA, 



1/2 







-M, 



1/2 



(2.6) 



This formulation will be preferred to study (|2.4p since the eigenvectors are now 

given as $^ = (^j, 0), $J = (Oj^j) where ^j are the eigenvectors of ^o, and 

A^f = ±iy/X]^f, with Ao'^j = Xj'ij- (At this step, remember that Aq is 

assumed to be a self-adjoint positive definite operator, so its spectrum is given 

by a sequence of positive real numbers going to infinity.) 

Now, setting 

-1 

T^Ao ] Ao, 



considering 



Ao.r 



I 



4/3- 1\ 2 



^fc+i/2 ^ ip';+^ ~ vl ^ ^^1/2 (^^+^^ 



fc+1/2 Lp';+^ - Vr 

y2,T = — - — 



^0,T 



,Al/^ ( ^1+^'^' 



(2.7) 



10 



system (|2.4p can be written as 

A. 



fc+l/2 fc-1/2 

Vr -yr 



fc-1/2 , fe+1/2' 



fee Z, 



with Ar = °^" 

I 



'iA. 



1/2 

0:1- 



(2.8) 



Under this form, one easily sees that another energy for solutions of (|2.4I) is 
given by 



2 



y^^+l/2 



X'^ 



"Pr 



fe+l 



fr 



Xo 



A, 



1/2 ( ^\ + y^ +^ 
fl,r 



(2.9) 



ATo 



As one easily checks, it turns out that the energies W'^^l'^ and E^^^l"^ (defined 
in (j2.5p ') are equivalent when working within a filtered class ^{5/t) at scale 1/r, 
thus making of no particular relevance to our purpose the fact that they do not 
coincide. 

Now, one can easily show that if y^/^ = $, where $ is an eigenvector of A 
given by (|2.6p corresponding to the eigenvalue ifi, then the solution j/t- of 
is given by 

y^+i/2 = exp (i^-^^kT I $, 



with /(a) = 2 arctan 



2 Vl + (/3 - l/4)a2 



, i? = oo. (2.10) 



Hence the Newmark approximation scheme (12.41) fits the assumptions of 
Theorem 11.31 

Also note that the observation ^11^9 + 621^9' for (|2.3p can be discretized in two 
different ways: 



A natural discretization consists in taking 



Bi 



Lp^ + (p'^+^ 



B. 



^k+i _ ^k 



for system (|2.4p . which corresponds to 



B „fc+i/2 ._ Ik, 4-1/2 ( k+1/2 _ k+1/2 



2^2 [Vi^r 



k+1/2 



k+1/2 
y2,r 



in the formulation 



A less natural discretization is as follows 



RA-'^xf:' (fi±!£i 



B2 



^k+i 



Vr 



which corresponds to 
By^^'l' := \b,A-,'I^ (fi:^^!^ - y\::^l^) + \b, [y^^ 
in the formulation 



/2 



fe+1/2 



V2,T 



(2.11) 



(2.12) 



(2.13) 



(2.14) 
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Note that Theorem 11.31 as stated, can only handle the second formulation 
(|2.13p in which the observation operator does not depend on r > 0, though it 
corresponds to a time-discrete observation operator of the form (|2.13p for (|2.4p , 
which seems less natural than (j2.1ip . 

Whether or not system (j2.4p is observable through (|2.1ip for general obser- 
vation operators Bi , B2 when the corresponding continuous system is observable 
is an open problem. 

Though, if 62 = 0, a trick allows us to get the same result as in Theorem ll.3l 
for an observation operator of the form 



Bi 



<y9^ + (^^+1 



Indeed, first apply Theorem II. 31 to the observation (|2.14p . There we obtain, for 
(5 > 0, r as in (|1.2ip . and initial data in €{S/t) back in the variable ipr- 



Vr-^T 



Xo 



A, 



1/2 (if°r+^l 



■0,T 



Xo 



<Cr Y. 

kT£{0,T) 



eu^"^Ai':; (^i^ 



Now, applying this identity to Af^ Aq / ipr, which is still a solution of (I2.4p 
with initial data in £((5/t). 



,1/2^-1/2/'^ 



^? 



Xo 



A ' 



l/2f^"r+vl 



Xo 



<Cr Y. 

kreiO.T) 



Bi 



VX + 'Pr+' 



But easy spectral computations show that, in the class £((5/r), there exists a 
constant C > depending only on b such that 



E^l"^ < C 



,1/2,-1/2 ( ^l-V'° 



Xo 



A ' 



l/2f'4+vl 



Xo, 



3 A representation formula, properties of the 
kernel pr and applications 

We will first recall basic facts on the discrete Fourier transform. We will then 
prove Theorem 11.11 and give some estimates on the kernel pr in (jl.Sp , which we 
use in Section 13.51 to prove Theorem 11.31 

3.1 Discrete Fourier transforms 

Let us introduce the definition of the discrete Fourier and inverse Fourier trans- 
forms: 
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Definition 3.1. Given any function Ut defined on rZ, we define its discrete 
Fourier transform J>[u^] at scale t as: 

J^^[M^](/iT-) = r\JuT(fcT)exp(— i^T-fcr), ^t- S (— tt/t, tt/t). (3.1) 

feez 

For any function v £ L^(— tt/t, tt/t), we define the inverse Fourier transform 
T^jT^lv] at scale r > as: 

J^^^[v]{kT) ^ —— / v{iJr)e^p{ifJ'TkT) diir, fc G Z. (3-2) 

According to Definition 13.11 one easily cliecks that these transforms are in- 
verse one from another, so that, 



J> "'"[J>[Mr]](fc''') = '"T-(fcr), fc e Z, 

•^r[^rM«]](Air) = v{flr), Mr G (-Tt/t, Tt/t), 



Similarly as for the continuous Fourier transform, we also have the following 
Parseval identity: 

^ r ' \^AUr]{f^r)\' d/ir - T ^ |l.. (fcr) ^ . (3.4) 

These properties will be used in the sequel. 

In the following, for a Hilbert space H, the space L^(rZ;_ff) is the set of 
discrete functions Ur defined on rZ with values in H endowed with the norm 

l|wr|li2(^z.H) = Tj^WMkr^H ■ 



k£l 



3.2 Proof of Theorem [TTT] 

Proof of Theorem \1.1\ Expand j/° as 



Then, for all k E Z, according to (J1.4D . 

J/r = E aj*jexp(i/ij>fcr), 

and thus y{t) defined by (|1.9[) can be written as 

y(i) = Ty^PritjkT) y. Gj^j exp{iiijrkT) 

fcez j, \t,,\<5/T 

= y. aj^jT2^ PT{t,kT) exp{i^j^rkT). 

j, \tij\<s/r fcez 

E «J^J-^-[P-(*'-)](-Mj,r). (3.5) 

i, IwI<5/t 
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According to (|1.8p . for all /^r G {^t^/ttTt/t), 

-^r[Pr(<, -jJlMr) = CXp X (AirTJ . 



Due to the definition of g, for all j such that |/ij | < (5/r, g{fij^rT) — fJ-jT- It then 
follows from (13.51) that 



y(i) = V aj<^.jexp{ifijt)xitJ-j^T 



T 



3, IA'jI<'5/t 

But the choice of x implies that 

j/(i) = ^ aj$j exp(i^ji). 

i, ImjI<'5/t 

Hence y{t) solves (jl.ip with initial data y". D 

Remark 3.2. When the parameter R € M^UJoo} in (|1.6p is smc/i iftai /(i?) = tt 
j/i? is finite or lini/ = tt, i.e. w/ie??, / is bijective from {—R, R) to (— tt, tt), i/iew 

oo 

there is no need of introducing a cut-off function to get Theorem \1.1[ To be 
more precise, we have the following result: 

Proposition 3.3. Under the assumptions of Theorem ll.li if we further assume 
that f is bijective from {—R,R) to (— 7r,7r), then we can take the function x in 
Theorem ] l.ll to be identically one. In other words, if for t > we define Prfiit, s) 
by 

PrAt, ^)^^ C^ exp (^ail^llF] e^t^^s ^^^^ (^^ ^) ^ ^2, (3.6) 



2tT J_^/r 

we have the following result: for all y^ £ €.{R/t) and y^ is the corresponding 
solution of (II. 3p . the function y(t) defined by (|1.9I) with p^fl instead of pr is the 
solution of (jl.ip with initial data y^ . 

Indeed, under the additional assumption that f is bijective from {—R, R) to 
(— tTjTt), the function g is defined on the whole interval (— 7r,7r). 

Remark that all the numerical schemes presented in Section \^ fit the as- 
sumptions of Proposition \3.3\ and thus Proposition \3.S\ applies for a wide range 
of numerical schemes. 

But even under this additional assumption, the localization properties of the 
kernel function prfi may be very rough and are not suitable to derive good es- 
timates on the time of uniform observability for solutions yr of (|1.3p as in 
Theorem \1.S\ 



3.3 Localization of p^ 

We now analyze the function pr in ([1 



Proposition 3.4. Let f be a smooth function describing the time discrete op- 
erator T-r as in (jl.4p . assume (|1.5p - (|1.6p - (jl.7p . and fix S £ (0,-R). 

Let £ > such that S + e < R, and choose the function x in Theorem \1.1\ 
supported in {—f{S + e), f{S + e)) and real-valued. 
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LetT >{) and {t, s) G (0, T) x R 6e such that 

t + e<s inf {/'(a)} or s sup {f{a)}<t-e. (3.7) 

\a\<S+e \a\<5+e 

Then for all n € N, there exists a constant Cn,e independent of (i, s) € (0, T) x R 
such that for all {t,s) satisfying (J3.7I) . 

n T-2n-l 

inf {|/'(a)s-t|}2» 

|ct|<<5+e 

where pr is the kernel Junction given by p.Sp . 
Proof. Recall that pr is given by (J1.8D . Hence 

PT{t,s) = - — / exp f - (a^s -5(a^)t) j x(ar)dar- (3.9) 

Remark then that 

- -^ ( ^^P ( 7 ("^'^ ~ 5(ar)i) 

— (s - 5'(ar)t)^ + ;^.g"(ar)t j exp i ^ {arS - g{ar)t) 

For (t, s) € (0,r) X R satisfying ((XT)) or, equivalently, 
t sup {5'(ar)} + e<s or s <t inf {5'(q;t)} - e, (3.10) 

|Q^|</(5+e) |ax|</(<5+e) 

for some e > 0, the right hand-side of this identity does not vanish, and then 
we can write, for all ar with \ar\ < f{S + e), 

exp ( - {arS - g{ar)t) j = -T^Gr{ar)-p^ ( exp ( - {urS - g{ar)t) 

"^^^^-^"-^^ (.s-.g'(a.)t)'+»r.g"K)f ^'-''^ 

Hence, using the fact that x is compactly supported in (— /((5 + e), /((5 + e)), we 
get 

Pr{'t-,s) = - — / -r^G'r(Q:r)-3-^ f exp f - (a^s-.g(a^)t) j j X(ar)rfar 

= - — / exp f - {arS - .9(a^)i) j — (Gr(ar)x(aT)) dar- 

= (-!)" ^^ / exp ( --(arS-g(ar)i) j ( ^(G^(aT)-) j x(aT)rfaT, 



where n G N and f -^{Grioir)-) ) denotes the operator ^^(Gr(aT)') iterated 
n times. 



15 



We finally remark that, since x is smooth and compactly supported on 
{—f{^ + ^)j/(<5 + £)) and due to the explicit form of Gr given by p. lip , for 
any n G N, there exists a constant C„ e such that for all ar and (i, s) satisfying 

(^(a(a.).))"xK) 



, *-^n,£ 



(s-5'(a,)t)2"- 
This immediately yields p.8p . D 

3.4 The transmutation operator 

We then prove that the transmutation operator is bounded as an operator from 

L2(tZ) in L^{R). 

Proposition 3.5. Let f be a smooth function describing the time discrete op- 
erator Tr as in (fL4| . assume (fT3)) -([lT6 l) - (fTT)) . and fix S £ (0, i?). 

Lei e > such that 6 + e < R, and choose the function x in Theorem \1.1\ 
supported in {—f{6 + s), f{d + e)) and real-valued. 

For T > 0, set X,- the transformation defined for discrete functions Wt com- 
pactly supported on rZ with values in some Hilbert space H by 

Tr{Wr){t) ^ T^ Pr{t,kT)Wr{kT), (3.12) 

fcgZ 

where p^ is the kernel function given by (jl.Sp . 

Then the operator T^- is bounded from L'^{t'L,H) to L'^{M.,H) and 

\\^rh(L2(rZ:HVL^(R.H)) < llxIL , / SUp /'(a). (3.13) 



\siiL^rZ;H);L^R;H)) < llxlloo./ ^Up /'(a) 

y |a|<<5+e 



Proof For Wr in L'^{tZ;H) and z G L'^{M.,H), both compactly supported in 
time, using Fubini's formulas, we compute 

{Ir{Wr)it),z{t))Hdt 

f (^T.T- r ' e-'3^'^--'>'/-x{l^rT)e'''-''^dprWr{kr),z{t)) dt 
^ r ' xii^rT) (r Vu;.(fcT)e^^-^^ / z{t)e'^^^^-^^'/^ dt) dp, 

X[Prr)(J^r[Wr\[-pT),z[ ) dp. 



2^, X{~l^TT)(Tr{Wr\{[lr),^A ^j) '^A't, 
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where z denotes the Fourier transform of z. Hence, usmg p.4p . 



[Tr{Wr){t),z{t))ndt 



< 



-^ j \\^T[Wr\{^ir)\\]l d^ir\ 



27r 



— tt/t 



\x{^^^rT)[' 



1/2 



aiP-rT) 



1/2 



d^lr 



H 



< ll^r|lL2(rZ:H) llxllooj SUp /'(«) ||^;|Il2(r.h) 
y |Q|<<5+e 



where the last hne is justified by the change of variable ^lt = g^fiT-r), which is 
valid due to Assumption (|1.7p and the fact that x is supported in the interval 

{-f{S + e)J{d + e)): 






difJ-rT) 



dfiT 



H 



(5+£)/r 



27r 



(S+e)/T 



\x{f{^^r))\'\\z{m„.n^^r)d^, 



(S+e)/T 



< 



< 



llxlloo ( sup f'{a)\— \\z{fi)\\H d^i 



(5+e)/r 



llxIlL ( sup /'(a))— / ||z(^)||^ d/i 



< 11x11 

This concludes the proof of p. 131) 



sup /'(a)) ||z||^2(R.ff). 

|q|<<5+6 ' 



U 



3.5 Proof of Theorem [TTHl 

Proof of Theorem \1.3[ Let y'^ e €{S/t) and y^ the corresponding solution of 
(|1.3p . Using Theorem ll.il and setting ?/(t) as in (II. 9p with x compactly sup- 
ported on {—f{6 + e),f{S + e)) for e > small enough so that S + e < R, we 
obtain the solution of p.ip with initial data y". 

Using the observability (|1.18p of the continuous system (|l.ip , we thus obtain 



To 



yll<C/ \\Bymtdt<C 



Tj2Pr{t,kT)By'^ 



feez 



dt. (3.14) 



Estimate (I3.14p is the center of our argument. Now, we only have to check that 
the right hand-side of (J3.14I) can be bounded by the right hand-side of (jl.20p . 
To do this, we set 



Ti,. 



To 



\ci\<o+e 



■ e and ti_. 



-e, 



(3.15) 
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and write 

Jo 
where 



Tj2Pr{t,kT)By^ 



kei 



dt < 30<i,,, + 30ti.,,Ti,. + 30>Ti,,, (3.16) 



O 



dt, 



ti.e.Ti,, = / ^ XI Pr{t,kT)By'^ 

° ti,,<fcr<Ti,^ 

"'0 I ,_-^^ 



di, 



fcr>Ti , 



(3.17) 
(3.18) 
(3.19) 



Now, using Proposition l3 .41 and ProDOsition l3.5l we prove the following facts: 
• There exists a constant C independent of r > such that 



O 



ti.s.Ti. 



< 



,ii2 



Gr E WrWu 



(3.20) 



ti,,<kT<Ti„ 



• For all n £ N*, there exists a constant C„ = Cn.e,p,5 independent of t > 
such that 

/O I/O ^ /^ 4n— 2— 2p |L,0||2 



lAT 



(3.21) 



Indeed, estimate (I3.20p can be deduced immediately from Proposition 13.51 
by choosing Wrikr) — By'; for fcr € {ti^^,Ti^^) and for fcr ^ (ii^e,ri^e). 
Let us then focus on p.2ip . According to Proposition 13.41 we have 



0<*,,. < / Ysup||%^-||^) fr E \Pr{t,kT)\] dt, 



(3.22) 



Using that the discrete semi-group (|1.3p is conservative in norms X and I?(^p), 
preserves £(5/t), and that S satisfies (|1.17p . we immediately have that 



a 



P'^ IL.0||2 



sup iSy^' < ||« II 



y 



(3.23) 



Besides, according to Proposition 13.41 for all t e (0,ro), and n > 1, 

r" T-2n-l 

r E IP^(^'fc-)l<- E 7 ^^""^ 

fcr<ti,^ fcr<ti,e ( inf |/'(q;)A:t — t| 



^^ E 



CnT 



2n-l 



krKtiAt- inf {f'(Q;)}A:T) 

V |a|<5+e / 



-<C„r^"~\ (3.24) 



for some constant C„ depending on n > 1. 



18 



Similar estimates can be done to bound C'>Ti.j , yielding p. 211) immediately. 
We now conclude the proof of Theorem II. 31 Estimates (I3.16P together with 
(|5:^ and ((X^ show that 



TY,Pr{t,kT)By'; 



kei 



,2 

\\u 



dt<CT J2 ll^y'l 

tl,,<fcT<Ti,e 

+ C„T4"-2-2j>||yO||^_ (325) 



Thus (l3TT4|) implies 



\y \\x 



tl,s<fcT<Ti,, 



.11 2 

u • 



(3.26) 



Since the discrete semigroup (|1.3p is conservative, we can shift the time in p.26p . 
and obtain 



\y'\\l{l-C^r^-'-'n<Cr ^ \\By';\ 



u ■ 



(3.27) 



0<fcT<Ti_e-tl,e 

Taking n > 1 + p and r > small enough, we obtain (|1.20p with 

To 



Te — Ti.e — tie — — r 



2e. 



This concludes the proof of Theorem 11.31 by taking e > small enough. D 



\a\<o+e 



4 A reverse formula and its applications 

4.1 A reverse formula 

In this section, our aim is to explain that the representation formula derived in 
Theorem 11.11 allowing to write the solutions of the continuous abstract equation 
(|l.ip as functions of those of the time-discrete one (II. 3p can be reversed: 

Theorem 4.1. Let f be a smooth function describing the time discrete operator 
Tr as in pri| . assume P3)) - ([T^ - pTT)) . and fix S £ {0,R). 

Let f be a smooth function describing the time discrete evolution as in (|1.4p . 
and assume (11.51) . p.6p and (11.71) . 

Let X be a C°° function compactly supported in {—R,R) and equal to 1 in 
(—6,5). For T > 0, define then qr(t,s) for {t,s) G R^ as 

'^^(*' '^^hi """P [^^^) X O^r) e-'^* dM. (4.1) 

Then, ify'^ G £((5/r) andy{i) is the corresponding solution of (jl.ll) . the function 
yr defined for k E "Z by 

y^= f qAt,kT)y{t)dt (4.2) 



is the solution of (jl.3p with initial data y^. 
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Proof. The proof is similar to the one of Theorem ll.il Write the solution y{t) 
of (|l.ip on the basis of eigenfunctions of A, and remark that for all /i G R 
satisfying |/u| < 5/t, due to the definition of qr in (|4.ip . for all s G K, 



/ 

Jr 



Qrit, s) exp(i/ii) at = XKP'''') cxp I = cxp ' 



Details are left to the reader. D 

Remark 4.2. Similarly as in Remark \3.Sl if R in (|1.6I) satisfies R — oo, 
Theorem. {jn\ still holds when taking % = 1, which consists in replacing q^- in 
(|4.ip by qT.o defined by 

<7t,o(^, s) = — / exp e f' d/z. 

^T^ Jm \ T J 

Note that this definition has to be understood in the sense of 2?'(R^) as the 
integrand is not integrable. 

Similarly, one can prove the following: 

Theorem 4.3. With the notations and assumptions of Theorem \4-.l\ ifO < 5i < 
62 < R andy^ G £(52/T)n£((5i/r)''", the same result holds for a smooth function 
X compactly supported in {—R,R) and equal to 1 in {—S2,—Si) U {Si,S2)- 

Remark 4.4. The notation _L in €^{52/t) n £((5i/t)^ stands for the orthogonal 
complement with respect to the X scalar product, so that £((52 /r) n €.{5i/t)^ 
stands for the set Span{^j : Si < fJ.jT < 62}- 



Proof. The proof of Theorem 14.31 is the same as the one of Theorem 14.11 D 

As we will see below, these transmutation formula also yield non-trivial 
informations, for instance when dealing with uniform hidden regularity results 
or the optimality of the time-estimate (|1.21l) in Theorem 11.31 

4.2 Uniform hidden regularity results 

In this section, we are interested in the admissibility property - also called hidden 
regularity property - for system (II. ip . 

To be more precise, if S G C{T>{A^)M) for some p G N and Hilbert space 
U, system (jl.ip is said to be admissible for B if there exist a constant Co and 
a time Tg > such that any solution y{t) of (II. ip with initial data y^ G T>{Ap) 
satisfies 

''\\By{t)\\l<Co\\y%. (4.3) 



Note that, when p = 0, i.e. B G C{X,U), this property is straightforward 
since the energy of solutions y of (jl.ll) is preserved. However, when p > 0, 
this property is not at all granted and comes from subtle properties of the 
system under consideration, requiring suitable assumptions on B and in partic- 
ular an adequate interaction of B with the free dynamics generated by A. The 
paradigmatic example of such situation concerns the wave equation in a bounded 
domain with homogenous Dirichlet boundary conditions observed through the 
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flux on the boundary. In this case indeed, the operator B is not bounded on X 
but is still admissible, see [IS]. This additional property then allows to define 
solutions in the sense of transposition for (|1.19p . see e.g. |49| . 

Let us also remark that by the semi-group property, it is straightforward to 
show that if the admissibility estimate (|4.3|) is true for some To with a constant 
Co, it is true for all T with constant C{T) = Co(l + lT/To\). 

Finally, we point out that the estimate (|4.3|) is the reverse of the observability 
estimate (I1.18|) . and we may therefore expect that the strategy developed for 
getting uniform observability estimates for time-discrete approximations of (|1.3p 
also applies in the context of admissibility. This is indeed the case: 

Theorem 4.5. Assume that B € Z{T>{A'^)M) for some p E N andU an Hilhert 
space, and that B satisfies (|1.17p with constant Cp. 

Assume that equation (|l.ip is admissible for B at time Tq with constant Cq, 
i.e. for all y" £ 'D{A^), the solution yit) of (|l.ip with initial data y^ satisfies 

Let f he a smooth function describing the time discrete operator T,- as in 
([L4)) . assume (fTSl) - (fO)) - ([Tj)) . and fix 5 e (0,i?). 

Then, for all time T > 0, there exists a constant C(T) so that for all t > 
small enough, solutions yr of (|1.3p lying in €{S/t) satisfy 

r E \Mu<c\\yll- (4.4) 

feTG(0,T) 

Similarly as in Theorem ll.31 Theorem 14.51 is derived by careful estimates on 
the kernel function q^ in Theorem 14. II 

Note however that Theorem 14.51 can be found in [TJ", Theorem 6.5]. There, 
it is proved using an equivalent form of the admissibility property (14.31) in terms 
of packets of eigenfunctions, in the spirit of [44] for the observability of waves 
(see also [49]). Though, we will provide a short proof of Theorem 14.51 bv using 
the kernel qr given by Theorem 14.11 to show the flexibility and efficiency of our 
strategy. The proof of Theorem 14.51 is postponed to the end of the section. 

We first show that the kernel function q^- given by Theorem 14. II is mainly lo- 
calized in some parts of M^ when wc choose the function x compactly supported 
in (-(5 - e,6 + e), for e e (0, i? - 5): 

Proposition 4.6. Let f be a smooth function describing the time discrete op- 
erator T^ as in (T^ . assume PTS)) - fT^ ~ P^ • and fix 5 £ {Q,R). Let e > 
.such that S + e < R, and choose the function x in Theorem \4.1\ supported in 

{-5-£,5 + £). 

LetT > and {t, s) G R+ x (0, T) be such that 

t>s sup {f'{a)} + e or t<s inf {f'{a)}~e. (4.5) 

\a\<S+e |Q|<<5+e 

Then for all n G N, there exists a constant Cn,e independent of (t, s) G M x (0, T) 
.such that for all {t, s) satisfying (j4.5p . 

(^ T-2TI-1 

\a\<d+e 

where qr is the kernel given in (14. ip . 
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Proof. Again, we only sketch the proof which can be done following the one of 
Proposition 13.41 

By a change of variable, similarly as in p.9p , rewrite qr as 

Qrit, s) = - — / exp ( - (/(a)s - at) ) x(a) da. 

As in p.lip . we then remark that, for (t, s) e M+ x (0,T) satisfying (I4.5P and 

|a|<(,5 + e)/T, 

exp f - (/(a)s -OLt)\ = -T^Fr(a)-^ \Qyj;}[- {f{a)s - at) 

1 
with Ft {a) 



{f'{a)s-t)^-iTf"ia)f 

The rest of the proof follows line to line the one of Proposition 13.41 D 

Proposition 13.51 also has a comiterpart: 

Proposition 4.7. Let f be a smooth function describing the time discrete op- 
erator T^ as in (iLil) . assume (IT3)) - (fOl) - ((TT7)) . and fix S e {0,R). Let e > 
such that d + e < R, and choose the function x in Theorem \4^.1\ supported in 
{-6-£,5 + e). 

For T > 0, set J^ the transformation defined for functions w{t) compactly 
supported on M with values in some Hilbert space H by 



J^{w){kT):^ I qr{t,kT)w{t)dt, (4.7) 

where Qr is the kernel given in (|4.ip . 



Then J^ is a bounded operator from L'^{M.;H) to L^(rZ;iJ); 

llyll 

\ciLHR,Hy,LHTl.;H)) - I . ^ ^,, ■ 

mt /'(a) 



(4.8) 



Again, the proof is similar to the one of Proposition 13.51 Details are left to 
the reader. 



Proof of Theorem \4.5\ Theorem |43] can be derived similarly as Theorem ll.3l bv 
applying Theorem 14. 11 Proposition 14.61 and Proposition 14.71 D 

4.3 Optimality of the time-estimate (ll.2ip 

Our goal here is to prove that the time-estimate (jl.2ip is sharp. 
Let us begin with the following result: 

Theorem 4.8. Let f be a smooth function describing the time discrete operator 
T.r as in ((TH), assume (fT3)) - ([T^ - (fTT)) . and fix 6 G (0,i?). 

Assume that there exist p > and a constant Cp > such that (I1.17P holds. 

Also assume that there exist a time Ti and a constant C such that for any 
T > 0, solutions yr of (|1.3p lying in <Z{5/t) satisfy (ll.20p . 
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Then, for any < Si < S2 < S, for any time T satisfying 

T>Ti sup {/'(a)}, (4.9) 

there exist positive constants C and tq > such that, for all r € (0, tq), all 
solutions y of (jl.ip with initial data y^ £ ([,{62 /t) D £((5i/t)^ satisfy (|1.18p . 

Proof. Let < (5i < (52 < (5, consider a smooth even function x compactly 
supported on {—62 — £, —^1 + e) U {Si — e,S2 + e) for some e £ (0, 5 — S2), and 
set Qr as in Theorem 14.31 

Then, similarly as in Proposition 14.61 one can prove that for {t, s) e M x 
(0,Ti) such that 

t>s sup {f'{a)}+e or i<s inf {/'(«)}-£, (4.10) 

(5i-e<|a|<(52+e (5i-e<|Q|<(52+e 

for all n £ N, there exists a constant C„ > such that 

C T^"^^ 
di — £<|Q|<(52+e 



Similarly, Proposition 14.71 still holds. 

Hence, following the proof of Theorem 11.31 and in particular of estimate 
(|3.16p . we obtain that 

2 /•Tisup,^_^<|„|<,^+j/'(a)}+e 

r Y. Wr\\u<C WBymldt 

fcrG(0,Ti) -'-" 

+ C„r4"'2-2p||yO||^_ (4_^2) 

Applying (|1.20p and taking r > small enough, (jl.lSp follows with T = 
Ti sup5^_^<|„|<5^+^. /'(a) + 2£. ^_^ 

Since e can be chosen arbitrary small, the observability estimate (|1.18p holds 
for any solution of (|l.ip with initial data in €.{S2/t) H £(^i/t)-'- and time T as 
in (US. D 



We are now in position to prove that the estimate (jl.2ip is sharp. In order 
to show this, we specify the abstract system (11.11) to the simplest case fitting 
the assumptions, namely the transport equation at velocity 1 on the 1-d circle 
denoted by S and identified with the interval (0, 1) with periodic boundary 
conditions. In this case, the equation reads: 

dty + d^y = 0, {t,x)eRxS. (4.13) 

This correspond to an operator A — —dx, defined on X ~ L'j,{§), the space of 
periodic functions of period one in L^{0, 1), with domain ^{A) = ^i(S)- 

We consider the observation operator By = 2/(0), which is continuous on 
T>{A) = H^{S) and takes value inU =R . 

Since the solutions of the transport equations (j4.13p can easily be solved with 
characteristics, we get y{t, x) = y'^{x — t). It is thus completely straightforward 
to show that 

[\y"{x)\^dx^ f\y{t,0)\'dt. (4.14) 

Jo Jo 
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In particular, applying Theorem 11.31 for some discretization scheme corre- 
sponding to / satisfying (|1.5p - (|1.6p - (|1.7|) . taking S < R, for ah T satisfying 

^ > • f A,/ u ' (4.15) 

mf {/'(a)} 

|q|<<5 

there exist a constant C and tq > such that for all r e (0, tq) and j/t solution 
of 

2/^-+i = exp(*/(-*^T))y^, fceZ, y? = y°, 

withyl = -a:r, X>(yl) =i?#(S), A' = i|(S), (4.16) 

with initial data y^ G €{S/t), 



kT£(Q,T) 



To prove the sharpness of the time estimate (|4.15p . we are thus going to 
show the following: 

Theorem 4.9. Assume that f satisfies (jl.5[) - p.6p - (|1.7p and take S < R. 
There is no time T > satisfying 

^ < ■ f s\'( M (4.18) 

mf {/'(a)} 

\a\<S 

and constants C > and tq > such that for all r € (0, tq), solutions y-r of 
(HT^ with initial data y° G €{S/t) satisfy (liTf)) . 

Proof. We argue by contradiction and assume that there exist a time T satisfy- 
ing (|4.18p and positive constants C > and tq > such that for all r G (0, tq), 
solutions j/t- of (I4.16P with initial data y^ G ^{S/t) satisfy (J4.17I) . 
Since T satisfies (j4.18p . we can find [Si, 52] C (0, S) such that 

T< ' 



sup f'{a) 



According to Theorem 14.81 choosing 



fel sup f{a), 1 , 

we get the existence of a constant C such that for all r > small enough, all 
solutions y of (|4.13p with initial data y° G £(^2/^) n €{Si/t)-'- satisfy 

1 /.f 

0/„m2 -i^ ^ ri \„.l4- nM2 



|y^(x)|^dx<C / |2/(i,0)|^dt. (4.19) 

"'0 

Now, we show that this cannot be true for T < 1. In order to do this, let us 
remark that the spectrum of the operator A — —dx with domain ^{A) = H], (§) 
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on <Y = L'i{E) simply is given by the Fourier basis ($j(a;) — exp(2ij'7rx))jgz 
and corresponds to the eigenvalues (i/ij = 2iJT:)j^z- 

Besides, using that solutions of (|4.13p are simply given by y{t, x) — y°{x~t), 



\y{t,0)\^dt 



\y°ix)\'dx, 



l-T 



so that (I4.19|) can be rewritten as 

nl 



\y°{x)\'dx<C \y°{x)\'dx. 

Ji-f 



(4.20) 



We thus have to prove that (I4.20p cannot be true uniformly with respect to 
T > for y° lying in €{S2/t) D (L{Si/t)^. 

This can be proved by an explicit construction as follows. We choose 6o G 
((^i: ^2), a smooth compactly supported function x with support in (—1, 1) and 
with unit i^(R)-norm, and xq G § such that xq G (0, 1 — T). 

We then set, for xq € S to be chosen later on, 



y°ix) = Vri/^x V^2j7r - -^ cxp(2zj7r( 



(x-xo)). 



(4.21) 



First, let us note that the coefficients of y^ in the basis ^j{x) — exp{2iJTTx) 
vanish for 

2j7r-- >^. (4.22) 



This implies in particular that for all r > 0, the sum in (j4.2ip is finite and thus 
makes sense, and that for r > small enough, y^ indeed belongs to €{62 /t) n 

e:(<5i/r)^. 

We can also compute the i^(S)-norm of y° by Parseval's formula: 



f\yl{x)\'dx = V^Yl 
■Jo jez 



X [ Vr^jn - 



-^,^ flxiarda^^, (4.23) 



since the sum is a Riemann sum. 

We now claim that y!^ is concentrated around xq. In order to show this, for 
a function v = v{j) defined for j G Z, we introduce the discrete Laplacian 



Adw(j) = v{j + 1) + v{j - 1) - 2v{j), 



and remark that 



-Adexp(2ijn{x — xq)) — exp(2zj7r(a; — a;o))4sin^(7r(a; — a;o)). 
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Thus, for e > and \x — xo\ € (e, 1 — e), 

(4sin (7r(x - xq)))" f7^ V V^/ ^ ^ 



isi 



(4sin^(^(.T-xo))) ^.g 



,i:.".,_^,,.(,(^,,,_i.)) 



g2ij7r(£c-a:o) 



SO that 

sup |2/°(:r)|<C„,eT"-l/^ (4.24) 

Ix — Xo|S(£,l — e) 



where we used that 

/ / _ AW 



and that the number of non- vanishing terms in the sum is of order r^^'^, see 

Km . 

In particular, considering (|4.23l) and (I4.24p with e > such that {x S 
§, ds{x, xo) > e} C (1 — T, 1), where ds is the geodesic distance on S, we have 
found a sequence j/" e €{62 /t) n €{5i/t)-^ such that 

/' \yl{x)\''dx<Cl,T^--''^ while f\y%x)\^ dx ^J- , 
thus contradicting (|4.20p . D 

5 Further comments 

5.1 Fully discrete approximation schemes 

Our approach also applies in the context of fully-discrete approximation schemes 
for (|l.ip . or more generally, to time-discrete approximations of a family of time 
continuous equations depending on a parameter. 
Following [13], we introduce the following class: 

Definition 5.1. Let p e N, Cp > 0, Co > and Tq > be constant parameters 
and define the set S{p,Cp,Co,To) as the set of elements {A,X ,B,U) such that: 

• X and hi are Hilbert spaces; 

• A is a skew- adjoint unbounded operator defined in X with dense domain 
"D^A) and compact resolvent; 

• S e £,{'D{A^),1J) and satisfies (|1.17p with constant Cp; 

• System (|1.1[) is observable through B in time Tq and satisfies (|1.18p with 
constant Cq- 

Theorem II. 31 can then be reformulated as follows: 
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Theorem 5.2. Let p eN, Cp > 0, Cq > and Tq > 0. Let also f be a smooth 
function describing the discretization process as in (jl.4p satisfying (|1.5p ~ (|1.6p - 
(|1.7p anrf /ix d G (0,i?). Then, for all T satisfying (jl.2ip . f/iere ea;ist positive 
constants C and tq > sitc/i t/iai /or an?/ r G (0, tq), solutions yr of (|1.3p Z^/m;; 
in<L(5/T) satisfy (ll.20p uniformly for {A, X,B,U) m 5(p, Cp, Co, Tq). 



In particular, Theorem l5.2l allows to decompose the study of the observabihty 
properties of a fuUy-discrete approximation scheme in two steps. 

Indeed, if {A,X,B,14) belongs to some S{p,Cp,Co,To) and X is an infinite 
dimensional vector space, the usual strategy to construct a fully discrete ap- 
proximation of y' = Ay is to first design a space semi-discrete approximation 
scheme. If ft, > denotes the space semi-discretization parameter, these approx- 
imations are defined on finite-dimensional Hilbert spaces X^ and can be written 
as y'f^ = Ahyu- Similarly, the observation, originally given by By on the time and 
space continuous model y' — Ay, is approximated by Buyh for some operator 
Bh defined on Xh and with values in some Hilbert space lAh approximating hi in 
some sense. 

Theorem 15.21 then states that, if we can find some p e N and constants 
Cp, Co, To such that for ah ft > 0, {Ah^Xh^BhMh) e 5(p, Cp, Co,ro) then, 
taking / satisfying the assumption (|1.5p - (|1.6p - (|1.7p and fixing 5 £ (0, i?), the 
solutions yr,h of the fully discrete schemes 

<r=exp(z/HAr))2/,V, A; e N, y^^vl 

satisfy 

kTe{0,T) 

provided y° G €h{S/T), where €h{S/T) is the vector space spanned by the eigen- 
functions of Ah corresponding to eigenvalues of modulus smaller than S/t. 

Thus, in many practical situations, proving observability properties for fully- 
discrete approximations of (jl.ip uniformly with respect to both space and time 
discretization parameters is reduced to proving observability properties for time 
continuous and space semi-discrete approximations of (jl.ll) . 

We do not give further details on this strategy as we have already developed 
it in [131 Section 5] - except for the estimate (ll.2ip on the time which is new. 

Also note that our strategy also applies to derive uniform admissibility es- 
timates for fully discrete approximations of (II. ip - where uniform means with 
respect to both space and time discretization parameters - provided uniform ad- 
missibility estimates are proved for the corresponding space semi-discrete and 
time continuous approximations of (jl.ll) - here, uniform means with respect to 
the space semi-discretization parameter. 

In both situations, our approach successfully reduces the study of observ- 
ability/admissibility issues for the time and space discrete approximations of 
(|l.ip to the study of the observability/admissibility properties for the underly- 
ing time-continuous and space semi-discrete approximations of (11.11) . for which 
a large literature is available, see e.g. (THl M M IMl I3S] and the survey articles 

[sniiis]- 
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5.2 Discrete Ingham inequalities 



As a by-product of our analysis, we proved the following discrete Ingham in- 
equalities: 

Theorem 5.3. Let I — N or TL and {fij)j^i be an increasing sequence of real 
numbers such that, for some 7 > 0, 



(5.1) 



inf l/ij+i -AijI >7- 

Let f be a smooth function describing the time-discrete operator T^- as in (|1.4p . 
assume pTS)) - PT^ - (fL7)l and fix S e (0, R). 
Then for all time 

T> ..%,, , (5.2) 

there exist two positive constants C and tq > such that for all r G (0, tq), for 
all {aj)j,=i e £'^{L) vanishing for j £ L such that \^j,j\t > S, 

2 



-y 



c 



a,\'<T 



jei 



E 

kTe(0,T) 






,if(ti.jr)k 



< 



Cj2\a,\ 



(5.3) 



Theorem l5.3l has to be compared with the results of [10], which derived a dis- 
crete Ingham lemma. Indeed, |40j gives assumptions on an increasing sequence 
{Xj)jei of real numbers under which one can guarantee for all r £ (0,to) and 






^r E 

fere(0,T) 



E"j 



iXj kr 



< 



^Ei 



(5.4) 



It is proven in [3D] that (15. 4p holds when assuming the existence of a gap 7 > 
such that 

inf|A,+i-A,|>7, (5.5) 

and provided that the sequence Xj satisfies, for some p £ (0, 1/2), 

2tt-tP 



sup|Afc - Xe\ < 

k.l T 



(5.6) 



When both (|5.5I) and (|5.6p hold, see [10], for any time T > 2tt/^, there exists a 
constant C > depending only on 7 and p such that (|5.4p holds. 

Of course, the result in [40j is thus very similar to Theorem 15.31 as one can 
check by setting 

Xj = iSlhll^ j (^ J such that |^j|r < S, 

which satisfies (j5.5p with 7 = 7inf|Q|<5- /'(a) and (|5.6p due to the assumption 
S <R. 

Nevertheless, our approach yields various generalizations of discrete Ingham 
inequalities from the ones known in the continuous setting, for instance the ones 
in [ii|29]. 

For completeness, we briefiy explain below how Theorem 15.31 follows from 
Theorem 11.31 and Theorem 14.51 
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Proof of Theorem \5.3[ Let us consider X = £'^{I), set $j = {aij)i(zj with aij — 
Si J, and define the operators A and B by A^j = ifij^j, and ,B<i>j = 1, which is 
continuous on 'D{A) since 1/(1 + ^^) is summable under the assumption (|5.ip . 
According to Ingham's lemma [20_, for all T > 27r/7, there exists a constant 
C > such that for all (oj) G i^{I), ' 






.f < r 

Jo 





■1 
dt < C ) ; |a 



2 



This can be rewritten as follows: for all y" — X^i %'J'j ^ ^k^)-i the solution y of 
y' — Ay with initial data y" satisfies 

^||yi;.(,)<^^|%(t)prfi<C||y°||;,(^^. 

Hence we can apply Theoreni ll.3l and Theorem 14.51 immediately to A and B and 
the corresponding discretizations described by (jl.4p and satisfying assumptions 
p.5p - (|1.6p - (|1.7p . Theorem 15.31 then follows from the explicit form of the solu- 
tions of y^+^ — Ti-y^ with initial data y° = X^i "j'^i' which is simply given by 

5.3 Weak observability estimates 

In this section, we briefly focus on the case of weak observability estimates. To 
be more precise, we consider an observation operator B E £(I?(.4p),Z^) and we 
assume the following: there exist a norm ||-||^, a time T^, > and a positive 
constant C^ such that for all solutions y of (11.11) with initial data y*^ G ^{A^), 

Vt<Cl r\\By{t)\\ldt. (5.7) 



This property is a weak observability property for the system (|l.ll) . Roughly 
speaking, this property appears as soon as solutions y of (jl.ip satisfy the fol- 
lowing unique continuation property: 

By(f)=Oon(0,T)^y = 0, (5.8) 

since then one can simply define ||-||^ as 

WvX^ f \\By{t)\\ldt. (5.9) 

Of course, defining ||-||^ as in (|5.9p does not provide any other information than 
the unique continuation property (|5.8p . 

When trying to quantify unique continuation properties, it is then natural 
to try to find a norm ||-||^ which can be compared to the norms constructed on, 
for instance, the iterated domains of the operator A. 

To be more precise, we introduce the family Xr of Hilbert spaces indexed by 
r G M as follows: for n G N, we set Xn = 'D{A"'); for r G R+, we take n G N such 
that r G [n,n -|- 1], and we define Xr as the interpolate between A'„ and (Yn+i 
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of order r — n. We then define AV for r < as the dual spaces of Xr {^o = '^ 
is identified with its dual). The corresponding norms ||-||^ on AV are then the 
following ones: for y = ^ Qj^ji 



1?^ 



3 



"J' 



We thus assume that there exist a constant r e R and Cr > 0, such that for 
all y E Xp, 

Il2/Il.<a||2/IL. (5.10) 

Note that this definition makes sense for all y € Xp since the weak observability 
estimate (|5.7p guarantees that ||y||^ is finite for y £ Xp. Actually, estimate (|5.7p 
and B € £,{'D{A^),L{) also imply r < p. Also note that when r > 0, inequality 
(|5.7p is stronger than the one in (|1.18l) . We are thus mainly interested in the 
case r < 0. 

We then have the following variant of Theorem 11.31 

Theorem 5.4. Assume that B E £{'D{Ap),U) for some p e N andU an Hilbert 
space, and that B satisfies ()1.17p with constant Cp. 

Assume that equation p.ip is observable through B at time T^ with constant 
Cw and norm \\-\\^, i.e. for all j/° € T^^A^), the solution yit) of (|l.ip with 
initial data y" satisfies (J5.7I) . Also assume that there exist a constant C^ > 
and r < p such that (|5.10p holds. 

Let f be a smooth function describing the time discrete operator Tt- as in 
((L4)) . assume (fTSl) - (fTB)) -([L7 | . and fix 6 e (0,i?). Then, for all 

T> . .f" ,, , (5.11) 

inf {/'(a)} 

\a\<d 

there exist positive constants C and tq > depending on f , T , 5, p, Cp, C^, 
Tjju, Cr, r such that for any r G (0,to), solutions yr of (|1.3p lying in €{S/t) 
satisfy 

hrWUcr E P*yr\\l- (5-12) 

fcre(0,T) 

Proof. The proof of Theorem 15.41 is very close to the one of Theorem 11.31 see 
Section [5751 Indeed, we first introduce e > 0, a function x compactly supported 
on {—f{S + e),f{S + e)) and taking value one on {—f{6),f{S)) and the kernel 
function pr given by Theorem 11.11 Then setting ti^^ — —e and Ti^^ = e + 
Tu,/inf|„l<5{/'(a)}, we get the counterpart of (I3.25P : for some constants C,Cn 
independent of r > 0, 



TypAt,kT)By^r 



k<£ 



dt<CT Yl W'^yrWu 



C„t4"-2-2p|I 01.2 



X ■ 



Using then Theorem ll.il and the fact that y{t) defined by (|1.9p is a solution 
of (jl.ip . we have 

l|y II* - '-'^ Z^ IP2/r||;^ +'-^n^ 11?^ lU' 

tl,e<fcT<Ti,j 
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instead of (j?:^ . 

One then uses that, since j/° g £((5/r), there exists a constant C > inde- 
pendent of r G (0, 1) and y° such that 



|2 
^ I f/T-' Il7;"l 



Oi,:^ / J Cr^'' ||yO||; if r <0, 



-^ " [ CWy^Wl if r>0. 

Using then (|5.10p . estimate (13.25^ imphes 

\\y°\\l (l - C„r4"-2-2P+2mi„{r,0}j < ^^ ^ ||gyfc||^ ^ (5^13) 

ti.E<feT<Ti.E 

Thus, taking n > I + p — min{r, 0} and r > small enough, we obtain (|5.12p 
with 

mf {/'(a)} 

|Q|<d+e 

This concludes the proof of Theorem 15.41 by taking e > small enough. D 

Among the typical cases fitting our assumptions, let us quote the case of the 
wave equation on (0, 1) observed from a point xq G (0, 1) \ Q. In that case, the 
equation reads: 

dttV - d.,^y = 0, for (t, x) G (0, T) x (0, 1), 

2/(t,0) = y(i,l)=0, forte(0,r), (5.14) 

(2/(0, •), 9*2/(0, •)) = (y^ yi) in L\0, 1) x i/-i(0, 1), 

and the observation is given by y(t,XQ). 

Equation (J5.14I) indeed fits the abstract setting of (jl.ip by setting 

with X = L^{0, 1) X H-\n), V{A) = i?o(0, 1) X l2(0, 1). 
and the point-wise observation operator is given for smooth Y by 

by^b(^] =yi(xo). 

Sobolev's embedding easily shows that B is continuous from 'D{A^^^^^) to M for 
any e > 0. 

Besides, expanding solutions y of (|5.14l) in Fourier, one easily checks 

y{t, x) — V 2 2_^ ( o,j exp(ij7rf) + bj exp(— zJTri) ) sin(j7ra;), 

where the coefficients (oj), (bj) can be characterized from the expansion of the 
initial datum {y^,y^): if {y^,y^) are given by 

y°{x) = V2^aj sm{JTTx), y ^ (x) = \/2 ^ /3j sin( jttx) , 
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the coefRcients {uj), {bj) are given by 

In particular, using Parseval's identity, one easily gets 

f \BY{t)\dt= f \y{t,xo)\'dt = 2^(|a,f + |fe,f)sin^(j^xo) 
"'o Jo jy^ 



J>1 ^ -^ 

Of course, if xq ^ Q, one easily checks that sin(j'7rxo) cannot vanish for j e N, 
hence unique continuation holds. In particular, this implies that the semi-norm 
defined for 

Y = f ^^ \ = ( v^Ej>iajSin(J7ra;) \ 

by 

j>i \ J / 

is a norm satisfying (j5.7p for T^, = 2 and C^, = 1. 

Besides, one easily checks that the norms \\-\\^ in that case are simply given 

by 

Thus, if one wants to apply Theorem 15 .41 one should verify condition (|5.10p . i.e. 
that there exist r > and C > such that for all j S N, 

(l + 0-^)2r/2<Csin20-^xo). (5.15) 

It turns out that condition (J5.15I) is satisfied for a large set of irrational numbers 
Xq ^ Q, that we will denote by S in the following. 

Indeed, S contains the irrational numbers xq G (0, 1) whose expansion 
[0, xi, X2, • • • ,Xn,---] as a continuous fraction is given by a bounded sequence 
{xn), see [231 P-23], for which the following property holds: there exists a positive 
constant C such that for all fc e N, 

mi{\qxo-p\}>-. (5.16) 

psz q 

Since condition (|5.16p is stronger than (j5.15p . we deduce in particular that S is 
uncountable. 

Besides, S also contains all the irrational algebraic numbers according to 
Liouville's theorem: If xq € (0,1) is an algebraic number of degree d on Q, 
there exists c > such that for all q G N, 

inf{|ga;o-p|} > -r. 
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Theorem 15.41 applies when xq € S, and yields for instance the following 
observability result, corresponding to the Newmark method (|2.4p with /3 = 1/4: 
for all ^ > and T > 2(1 + (5^/4), there exists a constant C > such that for 
all T > small enough, solutions yr of 



1 '„fe+i_o„fc^„fc-n_fl ' i (y^+i + 2y^ + y^- 



(y^i-2y^- + y^i)-a 



fc-i^ 



yr + 2/r 



„0 yl - Vr 



= 0, 

for {k,x) G Zx (0,1), 
y\ forxe (0,1), 



(5.17) 



with initial data 

i.y^^y^) £ (Span{\/2sin(j7ra;), with :;' G N satisfying jii < S/t} 
all satisfy 



2/0 
2/1 



<Cr ^ \y';ixo)\' 

fere(0,T) 



(5.18) 



We conclude this paragraph by pointing out that similar weak observation 
properties often appear in the models for wave propagation on networks ([5]). 
Nervertheless, condition (I5.10p is sometimes violated as in the case of the wave 
equation observed through a subset which does not satisfy the geometric control 
condition, see [13 [53] . 



6 Open problems 

6.1 Non- conservative time-discretization schemes 

In this article, we focused on time-discretization schemes that preserve the en- 
ergy of the solutions, which is a natural class since the continuous model (II. ip 
also preserves the energy of the solutions. 

However, in many situations, it may be interesting to consider dissipative nu- 
merical schemes, adding some numerical viscous effects to avoid instabilities. A 
very simple scheme introducing such dissipation properties is the Euler implicit 
method, which approximates (|l.ip as follows: 



,fc+i 



■Ay';^ 



fc > 0, 



y'r-y°- 



(6.1) 



If y = $ j , one easily checks that the solution y^ is given by 

k 



y'r 



1 



1 — irnj 



$, 



Thus, when considering only eigenvectors such that t\ii\ > So > 0, solutions y^ 
decay exponentially. In particular, this implies that the observability inequality 
(jl.20p may be not appropriate and should be replaced as follows: for all yr 
solutions of (|6.ip . 



IT/t\ 



X 



kreiO^T) 
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We emphasize that, though the observabihty inequahties (|1.20p and (|6.2I) are 
completely equivalent when the time-discretization schemes preserve the energy, 
it is not anymore the case when considering the Euler implicit method (j6.ip . 
and estimate (|6.2p is indeed weaker than the observability property (|1.20p in 
that case. 

This seems to indicate that a representation formula similar to (I1.9P would 
involve kernels similar to the ones obtained to link the wave equation to the 
heat equations, see the ones proposed in PH E] for instance. 

Let us rapidly explain what is the new difficulty occurring when considering 
non-conservative schemes by transposing the formal arguments we developed 
in the introduction to the Euler implicit equation (16. 1|) . Similarly as in (J1.12p . 
solutions 2/1- of (|6.ip with initial datum y^ as in (|l.lll) can be written as 

y^ = ^a,$,exp('*^^^fcr 

but this time / is a complex-valued function given by 

f[a) = arctan(a) + - log(l + a^). 

Hence, similarly as in (J1.13p . one can introduce the formal continuous version 
of j/r given by 

Zr{s) — 2_,'-''J^J "^^P ( * ' ^ 

3 ^ 

But this solution is only well-defined for s > 0, and we are thus led to look for 
some kernel function pr = p,-(i, s) such that for all /i G M and t Cz R, 

e'^'=l pr{t,s)exp(i^^^s) ds, (6.3) 

where the domain of integration has been modified into W^ instead of K as in 
(ll.lSp . Also note that, to derive observability properties for the time-discrete 
model (|6.ip . we actually need the identity (|6.3I) only for t in a bounded set of 
time and |/i| < 5/r, which might help deriving appropriate kernels. 

But we do not even know if equation (|6.3p is solvable. It may be the case 
but probably to the price of involving more singular kernels, as for instance 
in [T3] where the kernel was only defined for bounded sets of time, and allows 
to express solutions of the conservative model (wave-type model) in terms of 
solution of the dissipative models (heat-type model) . 

In some sense, this problem could also be thought as follows: given an 
observable conservative system, can we guarantee nice observability properties 
for viscous versions of it? With that respect, it is worth pointing out the works 
[71 [ini [2S] considering the controllability properties of the transport equation 

dty + d,y~ed,,y^O, (i,a:) e (0,r) x (0, 1), 
y{t,0)^v{t), y{t,l)^0, forte(0,r), ^"- ' 

with vanishing viscosity parameter e > 0, which illustrate the difficulties one 
encounters when considering control issues for vanishing viscosity systems (see 
also [5D] for another example). Indeed, to our knowledge, it is still not known 
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what is the best time T guaranteeing that the systems (|6.4p are uniformly (with 
respect to e > 0) null-controllable, though the critical time is expected to be 1, 
i.e. the time needed to control the underlying transport equation obtained by 
setting e = and dropping the boundary condition at x = 1. 

The complexity of this singular limit problem is one more evidence of the 
intrinsic complexity of passing from strongly dissipative dynamics to conserva- 
tive ones. But the reverse problem is simpler. Indeed, in [37] for instance (see 
also [H]), it was proved that the null-controllability and the observability of the 
heat equation 

yt-Ay = 0, 

can be obtained as a limit, when e — >■ 0, of the corresponding wave- like properties 
for the one-parameter family of wave equations 

£ytt + yt - Ay ^ 0. 

This is also in agreement with the results in f34| allowing to write down 
solutions of the heat equation in terms of solutions of the wave equation. In 
our context, this corresponds to writing the solutions of the time-discrete im- 
plicit Euler schemes in terms of the time-continuous conservative dynamics. As 
one can easily check, this can be done with the same formula as in Theorem 
14.11 This allows obtaining interesting results about the uniform (with respect 
of time-step) admissibility properties or the optimality of the possible uniform 
time-discrete observability results. But, unfortunately, the key issue of getting 
uniform observability results for the time-discrete dynamics out of the contin- 
ues ones, requires a transformation expressing the solution of the conservative 
dynamics in terms of the dissipative one, and thus requires further analysis. 

6.2 Variable time-steps 

In many applications, it is important to allow the time discretization parameter 
to change in an adaptive manner. The precise study of such case seems to be 
out of reach by using our method which strongly relies on the use of discrete 
Fourier analysis. 
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